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Abstract Use of high dynamic range (HDR) images and 
video in image processing and computer graphics appli¬ 
cations is rapidly gaining popularity. However, creating 
and displaying high resolution HDR content on CPUs 
is a time consuming task. Although some previous work 
focused on real-time tone mapping, implementation of 
a full HDR imaging (HDRI) pipeline on the GPU has 
not been detailed. In this article we aim to fill this gap 
by providing a detailed description of how the HDRI 
pipeline, from HDR image assembly to tone mapping, 
can be implemented exclusively on the GPU. We also 
explain the trade-offs that need to be made for improv¬ 
ing efficiency and show timing comparisons for CPU vs. 
GPU implementations of the HDRI pipeline. 


1 Introduction 

The use of high dynamic range imagery in computer 
graphics and image processing has gained popularity in 
recent years. This can be attributed to the increased re¬ 
alism and visual quality that is afforded by use of HDR 
data. Techniques such as image-based lighting, environ¬ 
ment mapping, and special effects such as realistic mo¬ 
tion blur and the well known bloom effect all produce 
improved results if they use HDR data instead of low 
dynamic range (LDR) data (Reinhard et al 2010). 

This increased demand for working with HDR con¬ 
tent is well matched with the capabilities of modern 
graphics cards. Currently all modern graphics hardware 
support floating point textures and renderbuffers. This 
allows programmers to directly feed in a floating point 
HDR image and process it on the GPU. 

It is, however, often the case that HDR images used in 
graphics applications are created offline using the CPU, 
or they are obtained as pre-made from external image 
databases (Debevec 1998). However, given the large num¬ 
ber of independent pixel operations required to create 
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an HDR image, the process of HDRI assembly is very 
suitable to be implemented on the GPU. Thus, the first 
goal of this paper is to demonstrate how to create an 
HDR image from a set of bracketed low dynamic range 
(e.g. JPEG) images directly on the GPU by using the 
OpenGL API. 

Due to the limitations of conventional display de¬ 
vices, it is not possible to display HDR imagery directly, 
although this may change in near future as HDR dis¬ 
plays enter the mainstream (Seetzen et al 2004; Akyiiz 
et al 2007). Instead, their dynamic range needs to be re¬ 
duced followed by quantization into an integer 8-bit per 
color channel data type before they can be shown on a 
display device. Algorithms that perform dynamic range 
reduction are called tone mapping (or tone reproduction) 
operators (TMOs), and they range from simple linear 
scaling to sophisticated multi-scale approaches that at¬ 
tempt to simulate the human vision (see Devlin (2002); 
Reinhard et al (2010); Banterle et al (2011) for excellent 
reviews). 

Similar to the HDR assembly process, most TMOs 
are comprised of a large number of independent pixel 
operations which render them suitable for a GPU im¬ 
plementation as well. One of the most popular TMOs 
that belongs to this category is the photographic tone 
reproduction operator (Reinhard et al 2002). Thus, the 
second goal of this paper is to demonstrate how both 
the global and local versions of this operator can be effi¬ 
ciently implemented by using OpenGL fragment shaders. 
Different from previous work, we will show that the im¬ 
plementation of this operator neither requires expensive 
convolution nor Fourier transform operations to compute 
local adaptation luminances. 


2 Related Work 

In this section, we review the previous work that deals 
with optimizing the HDR imaging pipeline. Cohen et al 
(2001) introduced the idea of HDR texture mapping on 
the GPU. As contemporary graphics cards at the time 
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Fig. 1: A bracketed sequence captured with a Canon EOS 550D/T2i digital camera. Each exposure is 1-fstop apart from the 
next exposure in the series. 


of the study did not support floating point textures, the 
authors proposed a technique to simulate HDR textures 
by using multiple 8-bit textures. Battiato et al (2003), 
on the other hand, provided a state-of-the-art report of 
the HDRI pipeline from HDR image creation to tone 
mapping. However, implementation of the pipeline on 
the GPU was not discussed. 

The idea of tone mapping on the GPU was intro¬ 
duced by several authors (Goodnight et al 2003; Artusi 
et al 2003; Goodnight et al 2005). In Goodnight et al 
(2003) and Goodnight et al (2005), the authors imple¬ 
mented Reinhard et al (2002)’s tone mapping operator 
using fragment shaders. To implement the local version 
of this operator, they have devised an efficient GPU 
based convolution operation. Furthermore, the authors 
have shown how to apply the method to time-varying 
sequences such as HDR videos. Artusi et al (2003), on 
the other hand, proposed a general framework to speed¬ 
up global tone mapping operators by effectively dividing 
the workload between the CPU and the GPU. 

A real-time tone mapping operator that also mod¬ 
els the perception effects was developed by Krawczyk 
et al (2005). In this work, the authors modeled several 
important effects such as visual acuity, glare, and lumi¬ 
nance adaptation. Later work implemented a Reinhard- 
like operator on FPGA architectures (Hassan and Car- 
letta 200Tb,a). 

To summarize, previous studies made significant con¬ 
tributions to achieve real-time performance in tone map¬ 
ping. In this work, however, we explain how the full HDR 
imaging pipeline, from image creation to display, can be 
implemented in real-time. Different from previous work, 
we also show how a local tone mapping operator that 
utilizes local adaptation luminances can be implemented 
without having to implement neither convolution nor 
Fourier transform based approaches on the GPU. 


3 Theory 

In this section, we will briefly explain the theory be¬ 
hind the HDR image generation and tone mapping. Their 



Fig. 2: The combined result of the sequence in Figure 1 into a 
single HDR image which is tone mapped using the technique 
described in this paper. 


GPU implementation will be discussed in the following 
section. 


3.1 HDR Image Assembly 

HDR images can be created in several ways: direct cap¬ 
ture, rendering, and multiple exposures technique are 
among the most commonly used ones. Direct capture 
may become the de facto way of creating HDR images 
in future, but currently it requires special hardware fur¬ 
nished with HDR sensors and therefore is not commonly 
used by most photographers. Furthermore, most such de¬ 
vices impose other restrictions such as limited resolution, 
long capture times, and lack of color support (Reinhard 
et al 2010). Rendering, on the other hand, is only suitable 
for computer generated HDR imagery. 

The multiple exposures technique allows photogra¬ 
phers to take a bracketed sequence of LDR images using 
a conventional digital camera, and then merge them into 
a single HDR image. Figure 1 depicts such a sequence of 
9 exposures with each exposure 1-fstop apart from the 
next one. In that each exposure is properly exposed for 
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a different region in the scene, the final HDR image con¬ 
tains details in both dark and light regions (Figure 2). 
Owing to the fact that this technique allows generation 
of HDR images using off-the-shelf cameras, it is a popu¬ 
lar choice among photographers. 

A single pixel, Jj, of an HDR image can be computed 
by using the following formula in the multiple exposures 
technique: 


N f 

h = Y. L 

i=l 




N 

X 

i= 1 


w(pij), 


( 1 ) 


where N is the number of LDR images, Pij is the value of 
pixel j in image i, / is the camera response function, w is 
a weighting function used to attenuate the contribution 
of poorly exposed pixels, and U is the exposure time of 
image i. One can obtain an HDR image by computing 
this equation for all pixels. 

In this equation, the inverse of the camera response 
function, / -1 , is used to linearize (i.e. degamma) the 
LDR images as they are typically captured in the non¬ 
linear sRGB color space. / -1 can be recovered directly 
from the bracketed sequence using response curve recov¬ 
ery algorithms (Debevec and Malik 1997; Mitsunaga and 
Nayar 1999; Robertson et al 2003), or it can be assumed 
to match the sRGB standard. We adopt the latter ap¬ 
proach in this paper to benefit from OpenGL’s sRGB 
texture support. 


3.2 Dynamic Range Reduction 

Standard display devices such as televisions and com¬ 
puter monitors are designed to display 8-bit per color 
channel integer input streams (although video cards that 
can output 10-bit and monitors that can display them 
have been in use for some time (AMD 2008)). Due to 
this limitation, HDR images and video cannot be di¬ 
rectly displayed on standard display devices. To display 
them, their dynamic range needs to be reduced followed 
by quantization into 8-bit integers. The algorithms that 
perform this task are called tone mapping (or tone re¬ 
production) operators 1 . 

To date, various tone mapping operators have been 
proposed each with a different approach to dynamic range 
reduction. TMOs are generally classified as global and lo¬ 
cal with global operators applying the same compressive 
function to each pixel while local operators changing the 
shape of this function (thus the degree of compression) 
based on the statistics of the local neighborhood around 
each pixel. 

One of the most popular TMOs that is commonly 
used in practice, and that ranks high in user studies, is 
Reinhard et al.’s photographic tone reproduction oper¬ 
ator (Reinhard et al 2002). This operator comes in two 
flavors, namely the global and the local operator. 

1 Quantization into 8-bits is not part of tone mapping, but 
it is a necessary step to create displayable images. 


3.2.1 Global Operator 

The global operator starts by computing the key of the 
scene which indicates its overall subjective brightness. 
The key is approximated by the log-average luminance 
(see Section 3.3 for color space conversions needed to 
obtain luminance from color and vice-versa), L w : 

L w = exp (-t; yy o g((5 + L w (x,y))y (2) 

x,y 


Here, L w (x,y) indicates the world luminance 2 of pixel 
(x, y) and S is a small offset added to avoid singularity 
that may occur at log(0) if black pixels are present in 
the image. The summation is performed across the entire 
image. 

Once the log-average luminance is computed, it is 
mapped to a user defined value, a, based on the desired 
subjective brightness of the scene. This is accomplished 
by: 


L(x,y) = ^-L w (x,y). (3) 

For most scenes illuminated by moderate lighting, a can 
be set to 0.18. To render darker scenes, it may be reduced 
to 0.09 or 0.045 (or less), and for lighter scenes it may 
be increased to 0.36 or 0.72 (or more). 

Once the image is scaled in this manner, the actual 
dynamic range compression is performed using a sig¬ 
moidal compression function: 


L d (x,y) 


Lp^y) 

1 + L(x,y)’ 


( 4 ) 


where Ld(x,y) represents the display luminance. While 
this equation is guaranteed to bring all pixels into a dis¬ 
playable range, some intentional burning in bright areas 
may be desired to create a more natural photographic 
look. The amount of burning can be controlled by a user 
defined parameter, L w hite : 


Ld(x,y) 


L(x,y)( 1 + fy 2j 'i) 

_ white 

1 + L(x,y) 


( 5 ) 


In this final equation, all luminance values greater than 
L w hite will be mapped to 1; that is they will burn out. 
If L w hi te is set to infinity, this equation will reduce to 
Equation 4. 


2 The subscript w indicates world luminance which may be 
in absolute or relative units depending on the calibration of 
the image. 
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3.2.2 Local Operator 


The local operator resembles the global operator in that 
tone mapping is performed via a similar formula: 


Ld(x,y ) 


L{x,y) 

1 + Vi(x,y,s)' 


(6) 


The difference, however, is that V\ represents the local 
adaptation luminance in the neighborhood around the 
pixel (x,y). The size of this neighborhood is controlled 
by the scale parameter, s. V can be computed as 


Vi(x, y, s ) = L(x, y) ® Ri(x, y, s ), 


where Ri is a Gaussian profile of the form 


Ri(x,y, s ) 


1 / x 2 + y 2 \ 

7T(ais) 2 6XP V ( aiS ) 2 J ' 


( 7 ) 

(8) 


To determine the appropriate scale, Reinhard et al 
(2002) propose to compute the difference of Gaussian 
convolutions at different scales, V\ and V^. When the 
difference between the two convolution results is above a 
threshold, the appropriate scale is found. This, in effect, 
computes the largest uniform region around each pixel, 
which serves as an adaptation region for that pixel. This 
can be formalized as: 


V{x,y,s) 


Vi(x,y,s) - V 2 (x,y,s) 
2 ^a/s 2 + Vi(x,y,s) 


( 9 ) 


where 0 is a sharpening parameter. Here the goal is to 
find the largest scale s m that satisfies: 


\V (x, t/, s m ) | < e, 


( 10 ) 


where e is a user parameter. Larger values give rise to 
larger adaptation neighborhoods. Reinhard et al (2002) 
suggests using 0 = 8.0 and e = 0.05 as default parame¬ 
ters. 

The photographic tone mapping operator poses two 
challenges for a GPU implementation. First, the log- 
average luminance of the whole image needs to be com¬ 
puted - an operation which is not GPU friendly. Sec¬ 
ond, local adaptation luminances need to be computed 
for the local operator. This amounts to convolving the 
image with filters of varying sizes, which is also not a 
GPU friendly operation. In this paper, we show that both 
problems can be solved by judicious use of mipmapping. 


3.3 Dealing with Color 

The dynamic range compression described in the previ¬ 
ous section expects luminance values as input. However, 
in practice, we typically deal with color images. To con¬ 
vert color values to luminance, we need to employ color 
space transformations. After tone mapping we can in¬ 
vert these transformations to retrieve the modified color 


values. In this section, we briefly highlight the key fea¬ 
tures of these color space transformations. For a more 
complete treatment, we refer the reader to literature on 
color imaging (Wyszecki and Stiles 2000; Reinhard et al 
2008). 

To compute the luminance value for a given color 
triplet, we first need to know its color space. If this in¬ 
formation is not available, we can assume that the HDR 
image is in the sRGB color space as this is the default 
output color space for most digital cameras. We also as¬ 
sume that the HDR image contains linear color values. 
This is also a reasonable assumption as the HDR gen¬ 
eration process typically linearizes the individual expo¬ 
sures before combining them into the HDR image. We 
can then convert an sRGB color value into its CIE XYZ 
representation with the following transformation (ITU 
2002 ): 


X' 


"0.4124 0.3576 0.1805" 


Rw 

Y w 

= 

0.2126 0.7152 0.0722 


G w 

z w 


0.0193 0.1192 0.9505 


B w 


In the CIE XYZ color space, the Y component encodes 
the luminance. Thus, Y w is equal to the world luminance 
L w that we used in the previous section. We can now 
compress Y w to obtain the display luminance Yd which 
is equal to Ld in Equations 4 and 5. 

The output RGB colors can be computed by: 

XiX (i2) 

where C = R, G, B and c is used for optional satura¬ 
tion adjustment. Setting c > 1 increases saturation while 
c < 1 decreases it. It is worth noting that all of these 
transformations described in this section are performed 
for each pixel independently, and thus are very amenable 
to benefit from GPU implementation. 


4 Mipmapping 

As mipmapping constitutes a key part of our algorithm 
which we use to compute the global average, L w , and 
local adaptation luminances, Vi and V 2 , a brief review of 
the concept can be useful. Mipmapping, first introduced 
by Williams (1983), is a commonly used technique to 
map texture images onto polygonal surfaces. The idea 
of mipmapping is to store a texture image as a pyra¬ 
mid of multiple levels, where each level contains a pro¬ 
gressively lower-resolution version of the original image. 
During texture mapping, the level which most closely 
matches the screen size of the polygon that is being tex¬ 
tured is chosen as the source image. 

In OpenGL, the mipmap levels for two dimensional 
textures can be explicitly provided by the programmer 
using glTexImage2D or glTexSubImage2D calls. Alter¬ 
natively, the programmer can request automatic gener¬ 
ation of mipmaps from the OpenGL server by using the 
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Level n Level n+1 



Fig. 3: A higher mipmap level is created from the lower 
mipmap level by filtered reduction. 



Mipmapping 


Fig. 4: The difference between proper convolution and 
mipmapping for computing local adaptation luminances. Ri 
indicate local adaptation regions at different scales. 


glGenerateMipmap function. In this case, each level of 
the mipmap chain is created from the previous level by 
using filtered reduction (the first level must be provided 
by the user). Although no specific filtering algorithm is 
enforced by the OpenGL standard, most implementa¬ 
tions use box filtering (Khronos Group 2012). Thus, each 
pixel in a higher mipmap level represents the local aver¬ 
age of pixels in the lower level (see Figure 3). 

A mipmapped 2D texture can be sampled in the frag¬ 
ment shader using the GLSL construct texture (s, xy, 
b). Here, s is a handler to the texture that will be sam¬ 
pled. xy indicates the coordinates inside the texture im¬ 
age, and b is a bias that will be added to the mipmap 
level computed by OpenGL. If the texture size is equal 
to the screen size of the polygon that is being rendered, 
one can think of b equal to the mipmap level index. 

As mentioned above, we use mipmapping to efficiently 
compute a measure of local adaptation luminance around 
each pixel. In the original algorithm of Reinhard et al 
(2002), this is performed by computing a Gaussian con¬ 
volution around each pixel. It is therefore appropriate to 
discuss the differences between these two approaches. In 
convolution, each pixel is placed in the center of a con¬ 
volution kernel and a local average is computed within 
that kernel. The kernel size can be increased to compute 
convolution over a larger neighborhood. In mipmapping, 
however, the downsampled versions of the original im¬ 
age are computed once using filtered reduction. Although 
this is very efficient as each pixel is used only once, it may 
give rise to an asymmetrical neighborhood for comput¬ 
ing local adaptation luminances. The difference between 
the two approaches is shown in Figure 4. 

In this figure, i?i, R 2 , and R% indicate the local adap¬ 
tation regions around the pixel shown in red. In Gaussian 
convolution, this pixel is always placed in the center of 
the convolution kernel. In mipmapping, this is not always 
the case as shown in the figure. We show in Section 6 that 
this difference has only a minor effect on the quality of 
the results, and therefore the heavy computational cost 
of convolution can be avoided in most cases. 


5 Practice 

In this section, we will demonstrate how the theory de¬ 
scribed in the previous section can be put into practice by 
using OpenGL. As it would be impractical to illustrate 
the entire implementation, we will focus on its most cru¬ 
cial features. In our implementation, we used OpenGL 
4.2 which was the latest version of OpenGL at the time 
of this writing. However, lower versions of the language 
can also be used as long as they support the required 
functionality such as mipmapping, floating point tex¬ 
tures, sRGB textures and framebuffers, and GLSL. All of 
these versions are supported in OpenGL version 2.1 with 
appropriate extensions and natively on 3.0 onwards. 


5.1 OpenGL Setup for HDRI Assembly 

To create an HDR image on the GPU, we need to access 
the pixels of the bracketed LDR images in the fragment 
shader. The most convenient way to achieve this is to 
upload LDR images as textures and sample from them 
using an appropriate sampler. The code snippet below 
demonstrates how to create these textures and the sam¬ 
pler: 

int numlmages = 10; // number of LDR images 
GLuint *ldrTex = new GLuint [numlmages]; 
GLuint IdrSampler ; 

glGenSamplers (1 , &ldrSampler ) ; 
glGenTextures (numlmages , IdrTex) ; 

glSamplerParameteri (IdrSampler , 

GL_TEXTUREJMINJFILTER, GL^EAREST) ; 

Listing 1: Generation of source sampler and textures. 

Note that we generate only one sampler as it will be 
shared by all of the texture units. Once the textures 
are generated, we can upload the LDR images which are 
stored as one dimensional arrays with color channels in¬ 
terleaved as red, green, and blue: 
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for (int i — 0; i < numlmages; H—|-i ) { 

glBindTexture (GL_TEXTURE_2D, IdrTex [ i ]) ; 
g lTexImage2D (GL_TEXTURE_2D, 0, GL_SRGB8, 
w, h, 0, GL_RGB, GL JJNSIGNEDJ3YTE, 
ldrlmg [ i ]) ; 

} _ 

Listing 2: Uploading of LDR images into textures. 

Here, w and h denote the dimensions of the LDR images. 
It is important to note that the internal format of the 
textures is set to sRGB. This will allow us to retrieve 
the linearized color values when we sample from these 
textures in the fragment shader. In other words, sam¬ 
pling from an sRGB texture will approximate the result 
of / _1 (pij) in Equation 1. 

We can now set up the source texture and sampler 
bindings. First, we need to bind the LDR sampler into all 
of the texture units as we want to use the same sampler 
for all units. Second, we need to bind each LDR texture 
into a different texture unit to be able to access them 
simultaneously in the fragment shader. These settings 
can be achieved by: 

GLint IdrSamplerUnits [ 1 6] ; 
for (int i — 0; i < numlmages; ++i) { 
IdrSamplerUnits [ i ] — i; 
glBindSampler (i , IdrSampler); 

gl Act iveText ure (GL_TEXTURE0 + i); 
glBindTexture (GL_TEXTURE_2D, IdrTex [ i ]) ; 

} _ 

Listing 3: Sampler and texture bindings. 

Here, note that in addition to binding textures and sam¬ 
plers, we initialize an array called IdrSamplerUnits with 
sequential integers from 0 to numlmages —1. This array 
will later be used to specify which sampler will fetch data 
from which texture unit in the fragment shader. 

The settings above complete the source texture and 
sampler setup. We can now perform the destination setup 
which is necessary to store the resulting HDR pixel val¬ 
ues. To achieve this, we can create a floating point tex¬ 
ture and attach it to one of the color attachment points 
of a framebuffer object (FBO), and make that FBO the 
current render target as shown in Listing 4: 

GLuint hdrTex , fbo ; 

glGenTextures (1 , MidrTex) ; 
glBindTexture (GL_TEXTURE_2D, hdrTex) ; 
glTexImage2D (GL_TEXTURE_2D, 0, GLJtGBA32F, 
w, h, 0, GLJtGB A, GLJTOAT, NULL); 
glGenFramebuffers (1 , &fbo); 

glBindFramebuffer (GLJ1RAWTRAMEBUFFER, fbo) ; 
glFramebuffer Text ure (GLJ1RAWTRAVIEBUFFER, 
GL_COLOR_ATTACHMENTO, hdrTex, 0); 


Listing 4: Destination setup to write out resulting HDR 
pixel values. 

The last operation we need to perform before the render 
call is to update the two uniform variables that will be 
used in the fragment shader. To this end, we first need 
to obtain the locations of these uniform variables, bind 
the HDR creation program, and then upload the values: 

// Get uniform lo cations 

GLint samplerLoc = glGetUniformLocation ( 
prgHDRCreate , ” IdrSampler” ) ; 

GLint imagesLoc = glGetUniformLocation ( 
prgHDRCreate, ” numlmages” ) ; 

// Switch to the HDR creation program 
glUseProgram (prgHDRCreate) ; 

// Initialize uniform variables 
glUniform 1 iv ( samplerLoc , numlmages, 
IdrSamplerUnits) ; 

glUniformli (imagesLoc , numlmages) ; 

Listing 5: Updating of uniform variables. 

It is important to note the usage of IdrSamplerUnits 
which was initialized with sequential integers in List¬ 
ing 3. By writing its value into the uniform array sam¬ 
pler variable IdrSampler, we establish a contract that 
in the fragment shader IdrSampler [0] will sample from 
texture unit 0, IdrSampler [1] will sample from texture 
unit 1, and so on. 

At this point we have completed all the necessary 
OpenGL API setup for HDR assembly. We can start the 
process by setting the viewport size equal to the image 
resolution, and drawing a quad to touch all pixels: 

// Draw a w by h quad to touch all pixels 
glViewport (0 , 0, w, h) ; 

DrawQuad () ; 

Listing 6: Draw call. 

This will initiate the execution of vertex and fragment 
shaders whose details are provided in the following sec¬ 
tion. 


5.2 Shader Setup for HDR Assembly 

The vertex shader that we need for HDR assembly is a 
simple pass-through shader which updates the position 
and texture coordinate attributes of each vertex: 

#version 420 

// Input position and texture coordinates 
layout (location =0) in vec2 inPos ; 
layout (location =1) in vec2 inTexCoord; 
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out vec2 texCoord; 
void main() { 

gl.Position = vec4(inPos, O.Of, l.Of); 
texCoord = inTexCoord ; 

} _ 

Listing 7: Vertex shader for HDR assembly. 

Note that this vertex shader is not specific for HDR as¬ 
sembly. In fact, we will use the same shader for tone 
mapping. The heart of the HDR assembly process is im¬ 
plemented in the fragment shader shown in Listing 8: 

#version 420 
in vec2 texCoord ; 

uniform sampler2D ldrSampler [ 1 6] ; 
uniform int numlmages; 

void main() { 

const int refid = numlmages / 2; 

float weightSum — 0.0; 

vec4 hdr = vec4(0.0, 0.0, 0.0, 0.0); 

for (int 1=0; i < 16; i++) { 

if (i < numlmages) { 

vec3 ldr = texture (IdrSampler [ i ] , 
texCoord ) . rgb ; 
float lum = luminance (ldr ) ; 
float w — weight (lum) ; 

float exposure = pow(2.0, i — refid); 

hdr. rgb += (ldr / exposure) * w; 
weightSum += w; 

} 

} 

hdr . rgb /= (weightSum + le—6); 

hdr . a = log (luminance (hdr . rgb ) + le—6); 

fragColor = hdr ; 

} _ 

Listing 8: Fragment shader for HDR assembly. 

The main function above calls two other functions namely 
luminance and weight to compute the luminance of each 
pixel and its contribution to the corresponding HDR 
value. Because we assume that the LDR images are cap¬ 
tured in sRGB color space, the computation of luminance 
is based on the ITU-R BT.709 primaries (ITU 2002): 

float luminance (vec3 color) { 

return color.r * 0.2126 + color.g * 0.7152 
+ color.b * 0.0722; 

} 


As for the weighting function, we need to use a function 
which attenuates the contribution of over- and under¬ 
exposed pixels while emphasizing the effect of properly 
exposed pixels. Several weighting functions are proposed 
in literature. We choose the tent function proposed by De- 
bevec and Malik (1997) due to its simplicity: 

float weight (float val) { 

if (val <= 0.5) return val * 2.0; 
else return (1.0 — val) * 2.0; 

} _ 

Listing 10: Weighting function 

This function assigns the highest weight for the pixels in 
the middle of the input range, and linearly decreases it 
for lower and higher pixel values. 

We note that Listing 8 closely adheres to the HDR 
assembly equation shown in Equation 1. The main dif¬ 
ference is that we assume the LDR exposures to be sep¬ 
arated by 1-fstop apart. This allows us to compute the 
exposure ratios in the pixel shader directly (note the use 
of refid), instead of getting them from the application. 
A second difference is that, we let the OpenGL do the 
linearization of LDR images for us by specifying an in¬ 
ternal format of sRGB as shown in Listing 2. If more 
accuracy is desired, the precomputed actual camera re¬ 
sponse can be provided to the shader through a uniform 
array variable. 

Finally, it is important to note that we write out the 
logarithm of the luminance into the alpha channel of 
the HDR image. This will be useful to compute the log- 
average luminance via mipmapping as explained in the 
next section. 


5.3 OpenGL Setup for Tone Mapping 

Once the draw call in the previous section completes, 
the HDR image will be stored in the texture hdrTex. 
For tone mapping, we can bind this as a source texture 
and sample from it to access the HDR color values. We 
can then perform dynamic range compression, and write 
out the resulting compressed pixel values into an sRGB 
texture to obtain the final displayable image. First let 
us demonstrate the generation of the tone map output 
texture, and its binding to the target FBO: 

GLuint tmTex; 

glGenTextures (1 , VtmTex) ; 
glBindTexture (GL_TEXTURE_2D, tmTex) ; 
g lTexImage2D (GL_TEXTURE_2D, 0, GL.SRGB8, w, 
h, 0, GL_RGB, GL_UNSIGNED_BYTE, NULL); 
glFramebufferTexture( GL_DRAW.FRAMEBUFFER, 
GL.COLOR-ATTACHMENTO, tmTex, 0); 


Listing 9: Computation of luminance 


Listing 11: Generation of tone map output texture. 
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We can now create a sampler to sample from the HDR 
image in the fragment shader. The reason that we cannot 
use the LDR sampler that we already created is that we 
need the HDR sampler to have mipmapping enabled. By 
sampling from the highest mipmap level we can obtain 
the log-average luminance of the HDR image which is 
needed for tone mapping. 

GLuint hdrSampler; 

glGenSamplers (1 , MidrSampler ) ; 
glSamplerParameteri (hdrSampler , 

GL_TEXTURE JV1IN .FILTER, 

GLTINEARJV1IPMAP NEAREST) ; 
glBindSampler (0 , hdrSampler); 
glActiveTexture (GL_TEXTURE0) ; 
glBindTexture (GL_TEXTURE_2D, hdrTex) ; 

// Create a mipmap chain 
glGenerateMipmap (GL_TEXTURE_2D) ; 

Listing 12: Sampler and texture setup for tone mapping. 

Finally, we can update our uniform variables that will be 
used in the fragment shader and draw a quad to initiate 
tone mapping. 


float key = 0.18 f ; 
float Ywhite = 1 e6f ; 
float sat = 1.0 f ; 

GLint hdrSamplerLoc = glGetUniformLocation ( 
prgTonemap , ” hdrSampler” ) ; 

GLint keyLoc = glGetUniformLocation ( 
prgTonemap , ” key” ) ; 

GLint YwhiteLoc = glGetUniformLocation ( 
prgTonemap, ” Ywhite”); 

GLint satLoc = glGetUniformLocation ( 
prgTonemap, ”sat”); 

glUseProgram (prgTonemap) ; 

// Update the uniform v ariables 
glUniform 1 i ( hdrSamplerLoc , 0); 
glUniformlf (keyLoc , key); 
glUniformlf (YwhiteLoc , Ywhite); 
glUniformlf ( satLoc , sat); 

// Enable sRGB framebuffer output and draw 
glEnable (GLTRAMEBUFFER_SRGB) ; 
glViewport (0 , 0, w, h) ; 

DrawQuad () ; 

Listing 13: Sampler and texture setup for tone mapping. 

Here, key, Ywhite, and sat are user-defined parameters 
and can be modified to change the appearance of the tone 
mapping result as will be demonstrated in Section 6. 


5.4 Shader Setup for Tone Mapping 

The vertex shader that we use for tone mapping is iden¬ 
tical to the vertex shader for HDR assembly (see List¬ 
ing 7). The main work for tone mapping is performed 
inside the fragment shader as shown below: 

#version 420 

uniform float key; 
uniform float Ywhite; 
uniform float sat ; 

uniform sampler2D hdrSampler ; 

in vec2 texCoord ; 

layout (location =0) out vec4 fragColor ; 
void main() { 

vec3 hdr = text ure (hdrSampler , texCoord). 
rgb ; 

float logAvgLum = exp (text ure (hdrSampler , 
texCoord , 20.0) . a) ; 

fragColor . rgb = tonemap(hdr, logAvgLum); 

} _ 

Listing 14: Sampler and texture setup for tone mapping. 

The tone mapping routine closely follows the description 
in Section 3.2. First the linear sRGB values are converted 
to XYZ. Tone mapping is then performed to compress 
the luminance. Finally, the compressed luminance is used 
to obtain the displayable RGB values with an optional 
saturation adjustment: 

vec3 tonemap(vec3 RGB, float logAvgLum) { 
vec3 XYZ = RGB2XYZ(RGB) ; 

float Y — (key / logAvgLum) * XYZ.y; 
float Yd = (Y * (1.0 + Y / (Ywhite * 
Ywhite))) / (1.0 + Y) ; 

return pow(RGB / XYZ.y, 

vec3(sat, sat, sat)) * Yd; 

} _ 

Listing 15: Global tone mapping implementation. 

The implementation of the RGB2XYZ routine is straight¬ 
forward and omitted for brevity. 

The tone mapping implementation in Listing 15 per¬ 
forms global tone mapping. For some applications, it may 
be desirable to perform local tone mapping as it better 
preserves the visibility of details. Previous GPU-based 
approaches for local tone mapping implemented convo¬ 
lution operations on the GPU. Here, we demonstrate 
that reasonable results can be obtained by simply using 
OpenGL’s mipmapping ability in lieu of convolutions. 
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vec3 tonemap_lc (vec3 RGB, float logAvgLum) { 

float La; // local adaptation luminance 
float factor = key / logAvgLum; 
float epsilon = 0.05, phi = 8.0; 
float scale [7] — float [ 7] (1 , 2, 4, 8, 16, 
32, 64); 

for (int i=0; i < 7; ++i ) { 

float vl = exp (texture ( hdrSampler , 
texCoord , i).a) * factor; 

float v2 = exp (texture (hdrSampler , 
texCoord, i+l).a) * factor; 

if (abs(vl — v2) / ((key * pow(2, phi) / 
(scale[i] * scale[i])) + vl) > 
epsilon ) { 

La = vl; 

break; 

} 

else 

La = v2; 

} 

float Yd = Y / (1.0 + La) ; 

} _ 

Listing 16: Local tone mapping implementation. 

To approximate convolutions, we compute V\ and V 2 
from consecutive mip levels. For this approach to work 
it is important to set the minification parameter of the 
sampler as GL_LINEAR_MIPMAP_NEAREST as shown in List¬ 
ing 12. Note that this approach is not only faster than 
computing convolutions as was done in Goodnight et al 
(2003), but also much easier to implement. 

Once the HDR image is created and tone mapped, 
the results can be downloaded back to CPU using the 
glGetTexImage function of OpenGL. 


6 Results 

In this section, we demonstrate representative results 
that were obtained by using the algorithms described in 
this paper. We will first show the effect of changing the 
tone mapping parameters on the resulting images, and 
then demonstrate that our GPU implementation pro¬ 
duces similar results to two reference CPU implemen¬ 
tations. We will then illustrate the performance that can 
be gained by using our method and then compare it with 
a standard convolution based approach. 

Figure 5 depicts tone mapped versions of two HDR 
images that were created using 9 exposures captured 
with a Canon EOS550D/T2i digital SLR camera. On the 
left column, we demonstrate the effect of changing the 
key parameter of the tone mapping operator. As it can 



Fig. 5: The left column depicts the tone mapping results with 
increasing key value in each row (0.18, 0.36, and 0.72 from top 
to bottom). The right column, on the other hand, depicts tone 
mapping results with decreasing burn-out threshold (10 6 , 5, 
and 2 from top to bottom). 


be seen, increasing the key value results in progressively 
brighter images. On the right column, we demonstrate 
the effect of changing the burn-out threshold, or L white in 
Equation 5. As this threshold is reduced we can see more 
pixels getting clamped at the highest possible value. For 
instance, while the details outside the window is visible 
in the top image, this region burns-out in the bottom 
image. Thus, this parameter can be used to controllably 
burn bright regions in an image to create an artistic ef¬ 
fect. An automatic method to estimate reasonable values 
for these parameters is explained by Reinhard (2003). 

We also illustrate the influence of the saturation pa¬ 
rameter in Figure 6. We remind that saturation adjust¬ 
ment is not part of tone mapping, but can be applied 
as a post processing operation to create the desired fi¬ 
nal look. As can be seen from the figure, setting a low 
saturation parameter such as 0.5 yields a more grayscale 
result while a high saturation parameter such as 1.5 ex¬ 
aggerates the color saturation. 

The difference between the global and local opera¬ 
tors is depicted at the top row of Figure 7. As expected, 
the local operator better preserves the visibility of de¬ 
tails as can be seen in the close-ups on the right. At the 
bottom row of the same figure, we show the results gen¬ 
erated by using a reference CPU implementation 3 . As 
the figure shows, our results are very similar to the CPU 
implementation. In fact, the visibility of the details on 


3 We used pfstmo_reinhard02 operator from pfstmo pack¬ 
age (Mantiuk et al 2007). 
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Fig. 6: The effect of changing the saturation parameter. The image on the left has the saturation parameter set to 0.5 and 
the image on the right to 1.5. The center image has no post tone mapping saturation adjustment (i.e. parameter set to 1.0). 


the book appears to have been better preserved by our 
method. The difference could be attributed to using dif¬ 
ferent scale factors. Whereas the reference implementa¬ 
tion uses 1.6 as the ratio of two scales, we had to use 2.0 
due to mipmapping. 

Next, we compare our results with two reference CPU 
implementations using a qualitative metric (Figure 8). In 
this figure, the top row shows the global photographic 
tone mapping results obtained by our method as well as 
the implementation of the same method in the pfstmo 
package (pfstmo_reinhard02) and the original implemen¬ 
tation of Reinhard et al (2002). On the second row, we 
can see the visible differences as detected by the dynamic 
range independent visual quality assessment metric (Ay- 
din et al 2008). Here, the green color indicates the loss 
of contrast, blue indicates the amplification of contrast, 
and red indicates the reversal of contrast. We can see 
that the differences between the two CPU implementa¬ 
tions are minor 4 and similar to the differences between 
our result and Reinhard et al (2002)’a original implemen¬ 
tation. Visual inspection of a selected region confirms 
this similarity. At the bottom two rows, we show the 
same result but this time for the local operator. Again 
the differences between the two CPU methods and our 
GPU method are comparable. The close-ups show the 
enhanced details. 

Table 1: Performance comparison of creating and tone 
mapping an HDR image on the CPU versus GPU in 
frames per second. Timings do not include disk I/O and 
GPU texture upload times. 


Device 

HDR gen. 

Global TM 

Local TM 

CPU 1 

0.18 fps 

0.12 fps 

0.015 fps 

GPU* 

65 fps 

137 fps 

103 fps 


1 Intel Core i7 at 3.20 GHz 

2 Nvidia GeForce GTX 590 


4 Such differences between two implementations can be 
caused by different post processing operations after tone map¬ 
ping such as normalization, clamping, quantization which are 
not elaborated in the original paper but are nevertheless used 
in the codes. 


We show further set of results obtained by our method 
together with the reference implementation in Figure 9 
using well-known HDR images. As can be seen from the 
figure, our results are qualitatively similar to the refer¬ 
ence CPU implementation, but are obtained at a fraction 
of the time of the latter. 

We provide a run time comparison to illustrate the 
performance benefits of our method. Table 1 lists the 
results of such a comparison obtained by creating and 
tone mapping an 18 megapixels (MP) HDR image cre¬ 
ated from 9 exposures captured by a Canon EOS 550D 
(Figure 1). In this test we used a high end CPU and 
a GPU. As the results indicate, both creating and tone 
mapping an HDR image on the GPU yields immense 
performance benefits. HDR assembly, on average, yields 
2 — 3 orders of magnitude improvement, while tone map¬ 
ping yields 3 — 4 orders of magnitude. If disk I/O and 
GPU texture upload times are included in the timings, 
creating an HDR image takes about 13.8 seconds on the 
CPU whereas it takes only 4.4 seconds on the GPU. 

As for the memory consumption, the total GPU mem¬ 
ory in bytes required for storing LDR images is given by 
N x w x h x 3, where N is the number of exposures, and w 
and h are the dimensions of the images. The HDR image 
occupies w x h x 4 x 4 bytes of memory as it needs to 
be in 4 component per pixel floating point format. The 
full mipmap chain requires approximately 1.33 times this 
number. 

We conducted more experiments to understand which 
part of the algorithm takes the most GPU time. As accu¬ 
rate measurement of timings on the GPU is not straight¬ 
forward, we omitted individual parts of our algorithm 
to see its effect on the frame rate. The results are re¬ 
ported in Table 2. We can see that the majority of the 
time is spent on texture look-ups from the source expo¬ 
sures. This is followed by the time it takes to generate 
a mipmap chain which approximately reduces the frame 
rate by 23%. Luminance and weight computations have 
very small impact on the performance. 

Finally, we investigated how long a standard convo¬ 
lution operation on the GPU takes. As can be seen in 
Table 3, when the kernel size is 7 x 7 or greater, the 
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GPU global 


CPU global 


GPU local 


CPU local 


Fig. 7: Global (left) vs. local (middle) tone mapping. As can be seen in the close-ups, the local operator better preserves the 
visibility of details. The top row shows the results obtained by our GPU implementation, whereas the bottom row shows the 
results of a reference CPU implementation (Mantiuk et al 2007). 


Table 2: Performance effect of the individual parts on HDR 
generation on the GPU. 


Configuration 

Frame rate 

Full 

No mipmap generation 

No texture look-up 

No luminance computation 
No weight computation 

65 fps 

84 fps 

107 fps 

69 fps 

69 fps 


Table 3: Performance of convolving an 18 MP image with 
varying sized kernels on the GPU. 


Kernel size 

Frame rate 

3x3 

210 fps 

7x7 

78 fps 

11 x 11 

54 fps 

15 x 15 

41 fps 

19 x 19 

33 fps 


convolution alone takes more time than our mipmap op¬ 
timized implementation. Given that typically larger ker¬ 
nels would be required to compute local adaptation lu¬ 
minances, the performance of the convolution approach 
is likely to be even lower in practice. This indicates that 
our algorithm is not only simpler to implement, but also 
outperforms convolution without compromising quality. 

These results underline the importance of transition¬ 
ing to a full GPU pipeline for both creating and tone 
mapping high resolution HDR images. 


7 Conclusions 

With high resolution HDR images becoming more com¬ 
mon in image processing and computer graphics appli¬ 
cations, their rapid processing is gaining importance. In 
this paper, we have shown how one can achieve real-time 
performances by implementing the full HDRI pipeline 
on the GPU. We demonstrated the feasibility of the ap¬ 
proach as well as the improved performance that it af¬ 
fords. We emphasized the key features of the implemen¬ 
tation to facilitate its reproduction by other researchers 
and programmers. While the full HDRI pipeline may 
contain other operations such as the camera response 
recovery, image alignment, ghost removal, etc., the skele¬ 
tal implementation provided here can serve as a basis to 
implement these other functionality as well 5 . 
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